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Collective  atomic  scattering  and  motional  effects 
in  a  dense  coherent  medium 


S.L.  Bromley1,  B.  Zhu1,  M.  Bishof1/  X.  Zhang1'"1,  T.  Bothwell1,  J.  Schachenmayer1,  T.L.  Nicholson1,1",  R.  Kaiser2, 
S.F.  Yelin3,4,  M.D.  Lukin4,  A.M.  Rey1  &  J.  Ye1 


We  investigate  collective  emission  from  coherently  driven  ultracold  88Sr  atoms.  We  perform 
two  sets  of  experiments  using  a  strong  and  weak  transition  that  are  insensitive  and  sensitive, 
respectively,  to  atomic  motion  at  1  pK.  We  observe  highly  directional  forward  emission  with  a 
peak  intensity  that  is  enhanced,  for  the  strong  transition,  by  >103  compared  with  that  in  the 
transverse  direction.  This  is  accompanied  by  substantial  broadening  of  spectral  lines.  For  the 
weak  transition,  the  forward  enhancement  is  substantially  reduced  due  to  motion.  Meanwhile, 
a  density-dependent  frequency  shift  of  the  weak  transition  (~10%  of  the  natural  linewidth) 
is  observed.  In  contrast,  this  shift  is  suppressed  to  <1%  of  the  natural  linewidth  for  the 
strong  transition.  Along  the  transverse  direction,  we  observe  strong  polarization  dependences 
of  the  fluorescence  intensity  and  line  broadening  for  both  transitions.  The  measurements 
are  reproduced  with  a  theoretical  model  treating  the  atoms  as  coherent,  interacting 
radiating  dipoles. 
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Understanding  interactions  between  light  and  matter  in  a 
dense  atomic  medium  is  a  long-standing  problem  in 
physical  science1,2  since  the  seminal  work  of  Dicke3. 
In  addition  to  their  fundamental  importance  in  optical  physics, 
such  interactions  play  a  central  role  in  enabling  a  range  of  new 
quantum  technologies  including  optical  lattice  atomic  clocks4  and 
quantum  networks5. 

The  key  ingredient  in  a  dense  sample  is  dipole-dipole 
interactions  that  arise  from  the  exchange  of  virtual  photons 
with  dispersive  and  radiative  contributions,  and  their  relative 
magnitude  varies  between  the  near-field  and  far-field  regimes. 
The  dispersive  (real)  part  is  responsible  for  collective  level  shifts 
and  the  radiative  (imaginary)  part  gives  rise  to  line  broadening 
and  collective  superradiant  emission6-8.  Intense  theoretical 
efforts  have  been  undertaken  over  many  years,  to  treat  the 
complex  interplay  between  the  dispersive  and  radiative 
dynamics9-18.  However,  experimental  demonstrations  that 
provide  a  complete  picture  to  clarify  these  interactions  have 
been  elusive. 

Collective  level  shifts  and  line  broadening  arising  from  the  real 
and  imaginary  parts  of  dipole-dipole  interactions  have  recently 
been  observed  in  both  atomic19-23  and  condensed  matter24 


systems.  The  modification  of  radiative  decay  dynamics  at  low 
excitation  levels  has  also  been  observed  using  short  probe 
pulses25-28,  and  interaction  effects  were  manifested  in  coherent 
backscattering29,30.  Although  simple  models  of  incoherent 
radiation  transport  have  often  been  used  to  describe  light 
propagation  through  opaque  media31,32  and  radiation  trapping 
in  laser  cooling  of  dense  atomic  samples33,  coherent  effects 
arising  from  atom-atom  interactions,  which  are  necessary  to 
capture  correlated  many-body  quantum  behaviour  induced  by 
dipolar  exchange,  are  beginning  to  play  a  central  role.  For 
example,  the  dipole-dipole  interaction  is  responsible  for  the 
observed  dipolar  blockade  and  collective  excitations  in  Rydberg 
atoms34-41;  it  may  also  place  a  limit  to  the  accuracy  of  an  optical 
lattice  clock  and  will  require  non-trivial  lattice  geometries  to 
overcome  the  resulting  frequency  shift42.  Previous  theoretical 
efforts  have  already  shown  that  physical  conditions  such  as  finite 
sample  size,  sample  geometry  and  the  simultaneous  presence  of 
dispersive  and  radiative  parts  can  play  crucial  roles  in  atomic 
emission10-13,43-45. 

In  this  work  we  use  millions  of  Sr  atoms  in  optically  thick 
ensembles,  taking  advantage  of  the  unique  level  structure  of  Sr  to 
address  motional  effects,  to  study  these  radiative  and 
dispersive  parts  simultaneously.  We  demonstrate  that  a  single, 
self-consistent,  microscopic  theory  model  can  provide  a 
unifying  picture  for  the  majority  of  our  observations.  These 
understandings  can  help  underpin  emerging  applications  based 
on  many-body  quantum  science,  such  as  lattice-based  optical 
atomic  clocks4,46,47,  quantum  nonlinear  optics39,  quantum 
simulations48  and  atomic  ensemble-based  quantum  memories49. 


Results 

Experimental  setup.  Bosonic  alkaline -earth  atoms  with  zero 
nuclear  spin  have  simple  atomic  structure  compared  with  the 
more  complex  hyperfine  structure  present  in  typical  alkali 
metal  atoms  that  complicates  the  modelling  and  interpretation 
of  the  experimental  observations.  For  example,  88Sr  atoms  have 
both  a  strong  lS0—lPi  blue  transition  (2  =  461nm)  and  a 
spin-forbidden  weak  —  3 Pi  red  transition  (7  =  689  nm),  with  a 
strict  four-level  geometry  (Fig.  la).  When  the  atoms  are  cooled  to 
a  temperature  of  ~  1  pK,  Doppler  broadening  at  461  nm  is 
~  55  kHz,  which  is  almost  three  orders  of  magnitude  smaller 
than  the  blue  transition  natural  linewidth,  T  =  32MHz.  To  an 
excellent  approximation,  atomic  motion  is  negligible  for  atomic 


coherence  prepared  by  the  461 -nm  light.  To  the  contrary,  the  red 
transition  with  a  natural  linewidth  r  =  7.5  kHz  is  strongly 
affected  by  atomic  motion.  By  comparing  the  behaviours  of 
the  same  atomic  ensemble  observed  at  these  two  different 
wavelengths  (Fig.  lb),  we  can  thus  collect  clear  signatures  of 
motional  effects  on  coherent  scattering  and  dipolar  coupling50,51. 

We  use  the  experimental  scheme  shown  in  Fig.  la,  to  perform  a 
comprehensive  set  of  measurements  of  fluorescence  intensity 
emitted  by  a  dense  sample  of  88Sr  atoms.  The  sample  is  released 
from  the  trap  and  then  illuminated  with  a  weak  probe  laser.  We 
vary  the  atomic  density,  cloud  geometry,  observation  direction 
and  polarization  state  of  the  laser  field,  and  we  report  the  system 
characteristics  using  three  key  parameters  as  follows:  the  peak 
scattered  intensity,  the  linewidth  broadening  and  the  line  centre 
shift.  For  example,  along  the  forward  and  transverse  directions  we 
observe  different  values  of  intensity  and  linewidth  broadening,  as 
well  as  their  dependence  on  light  polarization  (see  Fig.  lc). 


Figure  1  |  The  experimental  scheme  and  concept,  (a)  We  weakly  excite 
the  strontium  atoms  with  a  linearly  polarized  probe  beam  and  measure  the 
fluorescence  with  two  detectors:  one  in  the  forward  direction,  x,  and  the 
other  almost  in  the  perpendicular  direction,  z.  We  probe  two  different  J  =  0 
to  J'  =  1  transitions.  The  first  transition  is  a  ’’So  —  1P1  blue  transition  with  a 
natural  linewidth  of  r  =  32  MHz  and  the  second  is  a  ]So  —  3P-|  red  transition 
with  r  =  7.5  kHz.  (b)  In  the  coherent  dipole  model,  photons  are  shared 
between  atoms.  When  the  Doppler  broadened  linewidth  becomes 
comparable  to  the  natural  linewidth,  dephasing  must  be  considered.  At  our 
~1|iK  temperatures  the  Doppler  broadening  is  ^40  kHz,  meaning 
motional  effects  are  important  only  for  the  red  transition,  (c)  The  three- 
dimensional  intensity  distribution  predicted  for  a  blue  probe  beam.  The 
coupled  dipole  model  predicts  a  strong  103  enhancement  of  the  forward 
intensity  compared  with  other  directions  and  a  finite  fluorescence  along  a 
direction  parallel  to  the  incident  polarization.  The  speckled  pattern  is  due  to 
randomly  positioned  atoms  and  can  be  removed  by  averaging  over  multiple 
atom  configurations. 
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We  also  observe  motional  effects  on  the  red  transition  in  contrast 
to  the  same  measurements  on  the  blue  transition. 

In  the  experiment,  up  to  20  million  88Sr  atoms  are  cooled  to 
~  1  pK  in  a  two-stage  magneto -optical  trap,  the  first  based  on  the 
blue  transition  and  the  second  on  the  red  transition.  The 
atomic  cloud  is  then  released  from  the  magneto -optical  trap  and 
allowed  to  expand  for  a  variable  time  of  flight  (TOF),  which 
allows  us  to  control  its  optical  depth  and  density.  They  are 
subsequently  illuminated  for  50(100)  ps  with  a  large-size  probe 
beam  resonant  with  the  blue  (red)  transition  (Fig.  2a).  The 
resulting  scattered  light  is  measured  with  two  detectors  far  away 
from  the  cloud  (see  Fig.  la).  One  detector  is  along  the  forward 
direction  x  (detector  DF)  and  the  other  along  the  transverse 
direction  z  (detector  DT,  offset  by  ~  10°).  For  a  short  TOF,  the 
atomic  cloud  is  anisotropic  and  has  an  approximately  Gaussian 
distribution  with  an  aspect  ratio  of  Rx:Ry\Rz  —  2:2:1,  where  R{x,y,z} 
are  the  root-mean-squared  radii.  We  define  OD  as  the  on 
resonance  optical  depth  of  the  cloud,  OD  —  .2,  where  R± 

depends  on  the  direction  of  observation  with  R±tr  =  Rx  =  Ry  and 
R±,F  =  (RzRy)l/2  for  the  transverse  and  forward  directions 
respectively,  N  is  the  atom  number  and  k  is  the  laser  wavevector 
for  the  atomic  transition  (see  Supplementary  Note  1). 


other  hand,  the  wider  angular  area  of  enhancement  is  attributed 
to  the  longer  wavelength  of  the  red  transition.  The  forward 
intensity  strongly  depends  on  the  atom  number.  In  Fig.  2b,  we 
present  measurements  of  the  forward  intensity  Ix  versus  the 
transverse  intensity  Iz  at  a  fixed  0  =  2mRad  for  different  atom 
numbers.  The  intensities  are  normalized  to  those  obtained  at  the 
peak  atom  number  as  used  in  Fig.  2a.  To  the  first-order 
approximation,  the  transverse  fluorescence  intensity  scales  line¬ 
arly  with  the  atom  number.  Hence,  the  forward  intensity  of  both 
the  blue  and  red  transitions  scales  approximately  with  the  atom 
number  squared. 

In  the  forward  direction,  we  have  also  investigated  the 
linewidth  broadening  of  the  blue  transition  as  a  function  of  the 
atomic  OD.  By  scanning  the  probe  frequency  across  resonance, 
we  extract  the  fluorescence  linewidth,  which  is  found  to  be 
determined  primarily  by  the  OD  of  the  cloud  (open  squares  in 
Fig.  2c).  For  the  range  of  0  <  OD  <  20,  the  lineshape  is  Lorentzian 
(see  insets);  however,  the  observed  lineshape  starts  to  flatten  at 
the  centre  for  OD  >  20.  We  have  also  varied  the  atom  number  by 
a  factor  of  four,  and  to  an  excellent  approximation  the  linewidth 
data  are  observed  to  collapse  to  the  same  curve  when  plotted  as  a 
function  of  OD  (open  triangles). 


Forward  observations.  The  coherent  effect  manifests  itself 
most  clearly  in  the  forward  direction  (Fig.  2).  To  separate  the 
forward  fluorescence  from  the  probe  beam,  we  focus  the  probe 
with  a  lens  ( Lx )  after  it  has  passed  through  the  atomic  cloud  and 
then  block  it  with  a  beam  stopping  blade,  which  can  be  translated 
perpendicular  to  the  probe  beam  (Fig.  2a  inset).  The  same  lens 
(Li)  also  collimates  the  atomic  fluorescence  so  that  it  can  be 
imaged  onto  DF.  The  position  of  the  beam  stopper  can  be  used  to 
vary  the  angular  range  of  collected  fluorescence,  characterized 
by  the  angle  ( 6 )  between  x  and  the  edge  of  the  beam  stopper 
(see  Methods).  Using  the  maximum  atom  number  available  in 
the  experiment,  the  measured  intensity  IXto(6)  is  normalized  to 
that  collected  at  0max  =  7.5  mRad.  Both  the  blue  (square)  and  red 
(triangle)  transition  results  are  displayed  in  Fig.  2a.  For  the 
blue  transition,  we  observe  a  1,000-fold  enhancement  of  the 
normalized  intensity  for  6  <  0.5  mRad.  The  enhancement  is  also 
present  for  the  red  transition,  but  it  is  reduced  by  nearly  two 
orders  of  magnitude  at  small  6  due  to  the  motional  effect.  On  the 


Transverse  observations.  For  independent  emitters,  the  forward 
fluorescence  should  have  no  dependence  on  the  probe  beam 
polarization;  however,  the  transverse  fluorescence  (along  z) 
should  be  highly  sensitive  to  the  probe  polarization  and  it  is  even 
classically  forbidden  if  the  probe  is  z  polarized.  However,  multiple 
scattering  processes  with  dipolar  interactions  can  completely 
modify  this  picture  by  redistributing  the  atomic  population  in  the 
three  excited  magnetic  states  and  thus  scrambling  the  polarization 
of  the  emitted  fluorescence.  Consequently,  even  for  a  z-polarized 
probe  there  should  be  a  finite  emission  along  z  (see  Fig.  lc),  with 
an  intensity  that  increases  with  increasing  OD.  Our  experimental 
investigation  of  the  fluorescence  properties  along  the  transverse 
direction  is  summarized  in  Fig.  3.  Under  the  same  OD  along  z, 
the  y-polarized  probe  beam  (square)  gives  rise  to  a  much  more 
broadened  lineshape  for  the  blue  transition  than  the  z-polarized 
probe  beam  does  (triangle),  as  shown  in  Fig.  3a.  Meanwhile,  the 
peak  intensity  ratio  of  fypoi/Izpoi  decreases  significantly  with 
an  increasing  OD,  indicating  the  rapidly  rising  fluorescence  with 


Figure  2  |  Forward  scattering,  (a)  Comparison  of  forward  scattering  intensity  versus  angle  using  a  red  and  blue  probe  beam.  We  use  the  setup  shown  in 
the  inset,  to  block  the  probe  beam.  After  interacting  with  the  atoms  the  probe  beam  is  focused  using  a  lens,  which  also  collimates  the  fluorescence  from  the 
atoms.  We  block  the  probe  beam  using  a  beam  stopper,  which  we  translate  perpendicular  to  the  probe  beam,  to  change  the  angular  range  of  fluorescence 
collected  by  the  detector,  characterized  by  the  angle  (0)  between  x  and  the  edge  of  the  beam  stopper  (see  Methods).  The  measured  intensity,  lX/ o(9),  for 
each  probe  beam  is  normalized  to  the  intensity  at  6max  =  7.5  mRad.  The  dephasing  caused  by  motion  reduces  the  forward  intensity  peak  for  the  red 
transition,  (b)  Comparison  of  intensity  in  the  forward  direction,  IXI  versus  intensity  in  the  transverse  direction,  lz.  Both  are  varied  by  changing  N.  All 
measurements  are  made  at  6  =  2  mRad  (arrow  in  a)  and  normalized  to  the  intensity,  lXi0,  for  the  atom  number  used  in  a.  (c)  Linewidth  broadening  in  the 
forward  direction  measured  by  scanning  the  blue  probe  beam  detuning,  A,  across  resonance.  Example  lineshapes  for  different  ODs  are  shown  in  the  inset. 
Two  different  atom  numbers  are  used,  N  =  1.7(2)  x  107  (blue  squares)  and  N/4  (cyan  triangles).  The  dashed  line  represents  r  for  reference.  All  solid  curves 
are  based  on  the  full  theory  of  coupled  dipoles  and  the  band  in  c  is  for  a  ±20%  atom  number  uncertainty.  All  error  bars  are  for  statistical  uncertainties. 
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a  b  c 


Figure  3  |  Transverse  scattering,  (a)  Linewidth  broadening  for  the  blue  transition  in  the  transverse  direction  for  y  polarization  (open  squares)  and  z 
polarization  (open  triangles),  (b)  Intensity  ratio,  lypo\/lZpou  of  y  polarization  to  z  polarization  measured  in  the  transverse  direction  when  a  blue  probe  beam  is 
used.  For  low  optical  depths  single  particle  scattering  is  dominant  and  for  single  particle  scattering  almost  zero  intensity  is  predicted  for  z-polarized 
fluorescence,  as  this  polarization  points  directly  into  the  detector,  (c)  Linewidth  broadening  for  the  red  transition  in  the  transverse  direction  for  y-polarized 
light,  showing  a  similar  trend  to  the  blue  transition.  This  transition  is  more  sensitive  to  magnetic  fields;  thus,  a  large  magnetic  field  is  applied  to  probe  only 
the  m  =  0  to  m'  =  0  transition.  All  solid  curves  are  based  on  the  full  theory  of  coupled  dipoles  and  the  band  in  a,  b  and  c  is  for  a  ±20%  atom  number 
uncertainty.  All  error  bars  are  for  statistical  uncertainties. 


a  z-polarized  probe  when  OD  increases  (Fig.  3b).  For  the  red 
transition,  the  existence  of  Doppler  broadening  requires  the 
lineshape  data  to  be  fitted  to  a  Voigt  profile.  With  the  Doppler 
linewidth  dD  fixed  from  the  thermal  velocity  measured  in  free 
expansion,  the  Voigt  profile  determines  the  line  centre  and  the 
Lorentzian  linewidth  with  the  Gaussian  linewidth  determined  by 
the  temperature.  Figure  3c  displays  the  Lorentzian  linewidth 
obtained  with  a  y-polarized  red  probe  showing  a  strong  increase 
of  the  linewidth  with  OD. 


Spectral  broadening  and  shift.  To  a  good  approximation,  the 
dependence  of  the  linewidth  on  OD  along  the  forward  and 
transverse  directions  (for  the  classically  allowed  y  polarization  in 
the  single  scattering  limit)  is  similar.  However,  owing  to  the 
anisotropic  aspect  ratio  of  the  cloud,  for  the  same  TOF,  the  OD  is 
lower  along  z  than  along  x.  This  is  responsible  for  the  smaller 
broadenings  measured  along  z  than  along  x.  The  classically 
forbidden  polarization  direction,  on  the  other  hand,  exhibits  a 
different  scaling  with  OD,  which  is  understandable  given  that  the 
emission  in  this  case  comes  only  from  multiple  scattering  events 
with  dipolar  interactions.  The  transverse  linewidth  broadening 
for  the  red  transition  is  similar  to  that  of  the  blue,  and  it  does  not 
depend  sensitively  on  motional  effects.  This  behaviour  is  in  stark 
contrast  to  another  important  observation:  the  shift  of  the 
transition  centre  frequency.  Figure  4  contrasts  the  linecentre 
frequency  shift  observed  for  1S0  —  1Pi  (square)  and  1S0  —  3 Pi 
(triangle,  with  original  data  reported  in  ref.  52  and  see 
Supplementary  Fig.  1).  The  blue  transition  frequency  shift  is 
consistent  with  zero  at  the  level  of  0.004T  using  an  atomic  density 
of  1012cm-3.  However,  the  measured  density  shift  for  the  red 
transition  (normalized  to  the  transition  linewidth)  is  more  than 
one  order  of  magnitude  larger.  This  density- related  frequency 
shift  significantly  exceeds  the  predicted  value  based  on  general 
S-matrix  calculations  of  s-wave  collisions52  (2.18xl0_  10 
Hz  cm3  if  the  unitary  limit  is  used). 


Theory  model.  Before  we  turn  to  a  microscopic  model  to  obtain  a 
full  and  consistent  understanding  of  all  these  related  experimental 
observations,  we  note  that  semiclassical  models53  treating  the 
atomic  cloud  as  a  continuous  medium  of  an  appropriate 
refractive  index  can  give  an  intuitive  explanation  of  the 
linewidth  broadening  in  the  forward  direction.  Classically,  an 
incoming  electric  field  is  attenuated  as  it  propagates  through  the 
medium  according  to  the  Beer-Lambert  law  and  the  forward 
fluorescence  intensity  is  determined  by  the  same  mechanism. 
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Figure  4  |  Frequency  shift.  Comparison  of  frequency  shift  normalized  to 
the  corresponding  natural  linewidth  for  the  blue  and  red  transitions.  The 
blue  frequency  shift  is  consistent  with  0-0.004  of  T  at  an  atomic  density  of 
1012cm-3.  The  red  shift,  on  the  other  hand,  shows  more  than  0.1T  at 
densities  up  to  0.7  x  1012cm-3. 


This  simple  semiclassical  model  recovers  the  linear  dependence  of 
the  forward  width  for  small  OD  and  predicts  a  nonlinear 
dependence  of  the  linewidth  for  large  OD  and  a  flattening  of  the 
line  centre.  However,  we  find  that  this  semiclassical  approach 
cannot  provide  explanations  for  most  aspects  of  the  experimental 
observations. 

The  full  microscopic  model  builds  on  a  set  of  coherently 
coupled  dipoles.  Here,  each  four-level  atom  is  treated  as  a  discrete 
radiating  dipole  located  at  a  frozen  position,  but  coupled  with 
retarded  dipole  radiation,  and  it  is  driven  with  a  weak  incident 
laser  beam.  The  atomic  ensemble  follows  the  Gaussian 
distribution  observed  in  the  experiment  with  the  appropriate 
aspect  ratio.  By  solving  the  master  equation  in  the  steady  state,  we 
find  that  the  coherence,  bj=  Tr[| g)  (ej\ p],  of  atom  j,  located  at  r;-  is 
modified  by  other  atoms  as18,54-60: 


J  2(A 


A,f  +  r  £ 


a',m  yy 


Gqc.oc7  (*)’  I’m) 

(w“-9  ' 


(i) 


Here,  |^)=1S0,  |ea)  corresponds  to  the  three  excited  states 
of  lPi  or  3Pi ,  with  ae{x,  y,  z }  representing  the  Cartesian 
states.  In  addition,  p  is  the  reduced  density  matrix  of  the 
atoms  and  Syy  is  the  Kronecker  Delta.  The  driving  laser’s  linear 
polarization  state  £  is  along  y  or  z,  with  wavevector  k  along  x, 
Rabi  frequency  and  detuned  by  Aa  from  the  | g)  — ►  \ea) 
transition.  The  function  Ga>a'( r)  accounts  for  the  retarded 
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pairwise  dipolar  interactions  and  is  given  by  18,48,56 
Ga,o/(r)=  —  jp  [(^a,a'  —  +  (^a,a'  —  3fa7v)(^  —  ^he 

fluorescence  intensity  J(rs)=(E^(rs)E^  \rs)),  detected 
at  position  rs,  can  be  determined17,18  as  a  function  of  bj. 


2  r>2 


I(r .) 


9hzr 

16k2fi2rl 


Ztvb:, 


(brrs)(b;.rs)]eM^-»), 


j,m 


(2) 


with  /r  is  the  atomic  transition  dipole  moment  and  ks=/crs. 


Discussion 

To  understand  the  forward  enhancement  we  first 
consider  non-interacting  atoms  under  the  zeroth  order 
approximation.  The  atomic  coherence  is  driven  only  by 
the  probe  field  that  imprints  its  phase  and  polarization  onto 

the  atoms:  ^a(°)=  lSaf^_er^2,  where  A=Aa.  The  corresponding 


 9 h2r2 


-(N  +  iV2e“|ks“ko|2jR i*)  has 


intensity,  1(0=^^  4(^+r2/4)  ■ 

Lorentzian  profile.  It  also  exhibits  an  N 2  scaling  and  an  enhanced 
forward  emission  lobe,  with  an  angular  width  given  by  the  ratio 
between  the  transition  wavelength  and  the  transverse  size  of  the 


sample  A l/(kR±^)-  The  forward  lobe  reflects  the  construc¬ 
tive  interference  of  the  coherently  emitted  radiation  stimulated  by 
the  laser.  Outside  the  coherent  lobe  the  constructive  interference 
is  quickly  reduced  due  to  the  random  position  of  atoms28,59,61. 
The  longer  wavelength  of  the  red  transition  corresponds 
to  a  wider  angular  width  of  the  forward  lobe  for  the  red 
fluorescence. 

Simple  considerations  can  also  give  rise  to  a  qualitative 
understanding  of  atomic  motion-related  effects  on  forward 
enhancement.  Again  for  the  red  transition,  the  Doppler  effect 
introduces  random  phases  accumulated  by  5(j)~kv/r.  Here,  v  is 
the  thermal  velocity.  The  dephasing  reduces  coherent  photon 
emission  and  gives  rise  to  a  net  suppression  of  the  forward 
emission  intensity.  The  suppression  becomes  stronger  with  dDAT, 

with  ^d=a/8  ** qn2  the  Doppler  width.  Such  a  suppression  is 


clearly  observed  for  the  red  transition. 

To  address  the  linewidth  broadening  we  now  consider  atoms 
coupled  by  dipolar  interactions,  which  tend  to  emit  collectively  in 
an  optically  dense  cloud.  The  collective  emission  manifests  itself 
with  a  broader  fluorescence  linewidth.  Moving  to  the  first- 
order  approximation,  we  note  that  the  atomic  coherence 
acquires  contributions  not  only  from  the  probe  beam  but 

(0)  +  b-(1).  Here, 

0 


r 

'V 


also  from  the  surrounding  atoms,  with  bj  ^  ^ 

h?l)=J0kK^kX>  and <«'=  G“’“' (ri“  r™)eH' 

For  a  relatively  dilute  cloud  with  average  interparticle  distance 
f  1  /k,  the  far-field  interactions  dominate;  thus,  higher-order 
terms  beyond  1/r  can  be  neglected.  Dipolar  interactions  modify 
the  fluorescence  lineshape,  with  consequences  of  both  a  frequency 
shift  that  depends  on  the  cloud  peak  density  n0  and  a  line 
broadening  that  is  proportional  to  OD:  A  — >  A  +  A  and 


r  —  r  +  r,  with  A=-  3v/y  -r  and  r= f®  r.  Thus,  the 

first- order  approximation  provides  an  intuitive  picture  about  the 
role  of  dipolar  effects  on  the  lineshape. 

However,  in  a  cloud  with  an  increasingly  large  OD,  dipolar 
interactions  are  stronger  and  multiple  scattering  processes 
become  relevant.  The  first-order  perturbative  analysis  then  breaks 
down62-64.  The  full  solution  of  equation  (1)  based  on  the 
coherent  coupled  dipole  model  becomes  necessary  to  account  for 
multiple  scattering  processes  (see  Methods).  The  first  signatures 


arise  from  the  forward  fluorescence  intensity,  where  its  naive 
N 2  scaling  is  reduced  with  an  increasing  atom  number  as  a 
consequence  of  multiple  scattering  processes.  This  effect  is 
observed  in  both  red  and  blue  calculations,  and  it  is  expected  to 
be  more  pronounced  on  the  red  transition  due  to  its  longer 
wavelength.  However,  atomic  motion  leads  to  a  lower  effective 
OD,  which  tends  to  suppress  multiple  scattering  processes  and 
thus  helps  to  partially  recover  the  collective  enhancement. 
The  solid  lines  in  Fig.  2a, b  represent  such  quantitative 
theory  calculations  for  both  transitions,  which  agree  with  the 
experiment. 

Meanwhile,  for  the  linewidth  broadening  observed  in  the 
forward  direction,  it  becomes  evident  that  the  scaling  of  the 
linewidth  versus  OD  turns  nonlinear  at  large  values  of  OD.  The 
experimental  data  falls  within  the  shaded  area  in  Fig.  2c,  which 
represents  the  full  solution  with  a  20%  uncertainty  in  the 
experimental  atom  number.  Multiple  scattering  processes  are  also 
key  to  the  explanation  of  the  measured  fluorescence  along  the 
transverse  direction,  especially  for  the  classically  forbidden 
polarization  z.  Indeed,  for  both  intensity  and  linewidth 
broadening  observed  in  the  transverse  direction,  under  either  y 
or  z  probe  polarization,  the  full  model  (shown  as  shaded  areas 
in  both  Fig.  3a,b)  reproduces  well  the  experimental  results 
on  1S0  —  1Pi.  Taking  into  account  motional  dephasing 
(see  Supplementary  Note  2),  the  transverse  broadening  for 
—  3 P\  is  also  well  reproduced  as  shown  in  Fig.  3c. 

So  far,  we  have  shown  the  observed  effects  on  linewidth  and 
fluorescence  intensity  are  uniquely  determined  by  OD.  However, 
following  the  arguments  discussed  above,  the  frequency  shift 
arising  from  the  dipolar  coupling  is  expected  to  scale  with  atomic 
density,  |2i  |/T  oc  n0k~ 3,  which  includes  both  the  collective  Lamb 
shift  and  the  Lorentz-Lorenz  shift50,65.  For  our  experimental 
density,  this  effect  is  <10-3,  which  is  consistent  with  the 
observed  frequency  shift  for  the  blue  transition  (Fig.  4).  (It  is 
noteworthy  that  the  role  of  multiple  scattering  processes  is  to 
further  suppress  this  frequency  shift  mechanism50.)  In  contrast, 
for  the  red  transition  the  measured  density  shift  (normalized  to 
r)  is  significantly  larger  than  what  is  predicted  from  the  current 
treatment  of  interacting  dipoles;  it  is  also  much  bigger  than  the 
unitarity  limit  of  s-wave  scattering.  Qualitatively,  we  expect  that 
as  the  atoms  move  and  approach  each  other,  the  long-lived 
ground- excited  state  coherence  in  the  red  transition  can  be 
significantly  modified  by  the  collisional  process  and  open  higher 
partial  wave  channels.  We  can  thus  expect  a  larger  collisional 
phase  shift.  This  process  can  be  further  complicated  by  atomic 
recoil,  light  forces  and  Doppler  dephasing66. 

We  have  shown  that  a  coherent  dipole  model  describes  light 
scattering  in  a  dense  atomic  medium  with  collective  effects  and 
multiple  scatterings.  The  model  captures  the  quantitative  features 
of  the  experimental  observations.  Motional  effects,  as  manifested 
in  dephasing,  can  be  captured  in  the  model  as  well.  Our  results 
provide  useful  guides  for  further  developments  of  optical  atomic 
clocks  and  other  applications  involving  dense  atomic  ensembles. 


Methods 

Coherent  dipole  model.  Here  we  present  the  derivation  of  equation  (1).  The 
fundamental  assumption  is  to  treat  the  atoms  as  frozen  during  the  interrogation. 
This  is  an  excellent  approximation  if  hr  is  much  faster  than  other  energy  scales  in 
the  problem.  The  latter  condition  is  always  satisfied  in  the  case  of  the  blue 
transition.  For  the  /  =  0  to  /  =  1  configuration  exhibited  by  88Sr,  we  can  label  the 
/  =  0  ground  state  as  |g)  and  the  excited  /  =  1  states  using  a  Cartesian  basis 
10=10),  IO=(l  -  1)  -  I  +  l»/^2,  KH(I  -  1)  +  I  +  1»/VI  Here,  the  |0,  ±  1) 
states  are  the  standard  angular  momentum  ones.  In  the  Cartesian  basis, 
the  vector  transition  operator  for  the  j  atom  located  at  r;-  can  be  written  as 
bj=xbj  +ybj  +zbj.  Here  bJ=\g)j(ea\.  On  this  basis,  the  master  equation 
governing  the  evolution  of  the  reduced  density  matrix  of  the  N  atom  ensemble, 
p,  in  the  presence  of  an  incident  laser  beam  with  linear  polarization  £  can  be 
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s  =  - 1 E  Kt  [nibf  +  Qj  bi >  p]  - '  E  [s$  &? .  p] 


+  i 


J2A«[bfbJ,p]  +  (ibjpbft  -  {bfym,p}). 


(3) 


«,J8 


where  Q;  =Q^elk  r3  is  the  Rabi  frequecy  of  the  incident  field,  polarized  along  £ 

•  k=0)  and  detuned  by  Aa  from  the  atomic  transition  |g)  — >  |e“).  The 
parameters  g*£=gajj  (r;-  -  rm)  and  Jjjf  =/«,/*  (r;-  -  rm)  are  the  components  of  the 
elastic  and  inelastic  dipolar  interactions  between  a  pair  of  atoms  at  position  r y  and 
rm,  respectively,  and  are  given  by 


&,/i(r)  = 


3r 


To(^)  - 


kr 
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+  yiikr^rp 


+j2{kr)hrp 
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where  r=|rj  =  |r;-  —  rm|,  y„(x),  jn(x )  are  the  spherical  Bessel  functions  of  the 
second  and  first  kind,  respectively.  Here,  also  a,  j?  =  x,/orz  represent 
Cartesian  components.  The  symbol  <5yy  is  the  Kronecker  Delta.  In  the  low- 
intensity  limit,  we  can  project  the  density  matrix  into  a  state  space  including  the 
ground  state  |G)  =  \gi,g2,  •  •• gN )  and  states  with  only  one  excitation57-59  such 

as  (jot)  =  |^i,  ...e“,  ■  •  -gN ^ •  In  this  reduced  space,  the  relevant  equations  of  motion 
simplify  to 
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where  pG,c=  Tr[p|G}(G|],  pjaM=  Tr |p(fcf?E<“)]  and  pitt,G=  Tr[b“p], 

As  we  are  interested  in  the  situation  of  a  weak  probe  limit,  <C  T,  we  expand 
the  density  matrix  in  successive  orders  of  AT,  p=p^  +  +  p(2)  +  . . . ,  and 

keep  the  first-order  terms.  At  this  level  of  approximation,  Pg,g  =  C  Pja,mp  =  0  and 
only  the  optical  coherences  bj  =  f)p  G  evolve  in  time  accordingly  to  the  following 
set  of  linear  equations: 


dbj 

dt 


2QA 


E 


(10) 


Here,  Gay( r)  =  (fa,p(r)  +  ig^p(r))/r.  The  steady-state  solution  can  be 

dba 

obtained  by  setting  =0  and  then  solving  the  subsequent  3 N  linear  equations. 


Measure  the  enhancement  of  forward  fluorescence.  To  measure  the  scattered 
light  in  the  forward  direction,  we  use  the  setup  shown  in  the  inset  of  Fig.  2a,  to 
tightly  focus  and  block  the  probe  beam,  while  still  collecting  most  of  the  atomic 
fluorescence  on  the  CCD  (charge-coupled  device)  camera.  We  focus  the  probe 
beam,  after  it  interacts  with  the  atoms,  to  a  small  spot  with  15  pm  root-mean- 
squared  radius  and  block  it  using  a  beam  stopping  blade.  We  then  translate  the 
beam  stopper  perpendicular  to  the  probe  beam  by  a  distance  Ax  from  our  reference 
point  of  x  =  0,  which  we  define  as  the  position  of  the  beam  stopper  where  we  see 
the  greatest  fluorescence  without  saturating  the  CCD  camera  with  the  probe  beam. 
As  only  the  forward  direction  is  particularly  sensitive  to  positional  changes,  we 
convert  this  change  in  position  to  a  change  in  angle  simply  using  6=  arctan  yA^, 


where  the  first  lens  with  a  15-cm  focal  length  collimates  the  fluorescence.  In 
numerical  calculations,  the  CCD  camera  is  simulated  as  a  ring  area  centred  around 
the  forward  direction  and  the  average  intensity  collected  over  the  ring  is 
determined.  The  external  radius  is  set  to  be  large  enough  to  reach  the  angular 
region  outside  the  interference  cone  and  the  inner  angular  radius  0sim,  simulating 
the  blocking  of  the  signal  by  the  beam  stopper,  is  varied  accordingly  to  the 
experiment.  To  account  for  the  difference  between  <7sim  and  the  experiment  cloud 
size,  0sim  is  rescaled  so  that  we  satisfy  the  experimental  observation  that  at  0max  the 
enhancement  factor  drops  to  1. 
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